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charging energy. A 7r-transition of the supercurrent can be driven by tuning gate or 
bias voltages. 
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1. Introduction 

Non-equilibrium transport through superconducting systems attracted much interest 
since the demonstration of a Superconductor-Normal-Superconductor (SNS) transistor 
PQ. In such a device, supercurrent suppression and its sign reversal (7r-transition) 
are achieved by driving the quasi-particle distribution out of equilibrium by means of 
applied voltages [2j [3j HI [5] . Another interesting issue in mesoscopic physics is transport 
through quantum dots attached to superconducting leads. For DC transport through 
quantum dots coupled to a normal and a superconducting lead, subgap transport is due 
to Andreev reflection [61 [TJ El [9J QUI EH]. Also transport between two superconductors 
through a quantum dot has been studied extensively. The limit of a non-interacting 
dot has been investigated in [12] . Several authors considered the regime of weak tunnel 
coupling where the electrons forming a Cooper pair tunnel one by one via virtual states 
P31 El [15]. The Kondo regime was also addressed p31 QH OH HI UJ]. Multiple 
Andreev reflection through localized levels was investigated in [20], |2T] . Numerical 
approaches based on the non-crossing approximation [22] , the numerical renormalization 
group [23] and Monte Carlo [24] have also been used. The authors of [25] compare 
different approximation schemes, such as mean field and second-order perturbation in 
the Coulomb interaction. In double-dot systems the Josephson current has been shown 
to depend on the spin state of the double dot [26]. Experimentally, the supercurrent 
through a quantum dot has been measured through dots realized in carbon nanotubes 
[27] and in indium arsenide nanowires [28J. 

In this Letter we study the transport properties of a system composed of an 
interacting single-level quantum dot between two equilibrium superconductors where 
a third, normal lead is used to drive the dot out of equilibrium. A Josephson coupling 
in SNS heterostructures can be mediated by proximity-induced superconducting 
correlations in the normal region. In case of a single- level quantum dot, superconducting 
correlations are indicated by the correlator (di(0)d^(t)) , where d a is the annihilation 
operator of the dot level with spin a. To obtain a large pair amplitude, i.e. the equal- 
time correlator (d^), at least two conditions need to be fulfilled: (i) the states of an 
empty and a doubly-occupied dot should be nearly energetically degenerate and (ii) the 
overall probability of occupying the dot with an even number of electrons should be 
finite. For a non-interacting quantum dot, i.e. vanishing charging energy U for double 
occupancy, this can be achieved by tuning the level position e in resonance with the 
Fermi energy of the leads, e = [12]. In this case, the Josephson current can be viewed 
as transfers of Cooper pairs between dot and leads and the expression of the current 
starts in first order in the tunnel coupling strength T. 

The presence of a large charging energy U ^> kgT, T destroys this mechanism since 
the degeneracy condition 2e+U ~ is incompatible with a finite equilibrium probability 
to occupy the dot with an even number of electrons. Nevertheless, a Josephson current 
can be established by higher-order tunnelling processes (see, for example, [TB I [T4 l [15]). 
associated with a finite superconducting correlator (c^(O)df(t)) at different times. The 



N on- Equilibrium Josephson and Andreev Current through Interacting Quantum Dots 3 



amplitude of the Josephson coupling is, however, reduced by a factor T/A, i.e., the 
current starts only in second order in T, and the virtual generation of quasiparticles 
in the leads suppresses the Josephson current for large superconducting gaps A. In 
particular, it vanishes for A — > oo. 

The main purpose of the present paper is to propose a new mechanism that 
circumvents the above-stated hindrance to achieve a finite pair amplitude in an 
interacting quantum dot, and, thus, restores a Josephson current carried by first-order 
tunnel processes that survives in the limit A — > oo. For this aim, we attach a third, 
normal, lead to the dot that drives the latter out of equilibrium by applying a bias 
voltage, so that condition of occupying the dot with an even number of electrons is 
fulfilled even for 2e + U ^ 0. 

We relate the current flowing into the superconductors to the nonequilibrium 
Green's functions of the dot. In the limit of a large superconducting gap, A — > oo, 
the current is only related to the pair amplitude. The latter is calculated by means 
of a kinetic equation derived from a systematic perturbation expansion within real- 
time diagrammatic technique that is suitable for dealing with both strong Coulomb 
interaction and nonequilibrium at the same time. 

2. Model 

We consider a single-level quantum dot connected to two superconducting and one 
normal lead via tunnel junctions, see figure [TJ The total Hamiltonian is given by H — 




Figure 1. Setup: a single- level quantum dot is connected by tunnel junctions to one 
normal and two superconducting leads with tunnelling rates Tn and Tg L R , respectively. 

H D + J2 v =n,s l ,s r (H v + H tunnjV ). The quantum dot is described by the Anderson model 
Hd = J2a td\d a + Un^rii, where n a = d^d^ is the number operator for spin a =|, |, e is 
the energy level, and U is the charging energy for double occupation. The leads, labeled 
by T) = N, S L , S R , are modeled by H v = J2ka ^kc\ ka c vk(J - £ fc (A^cj^cJ^ + H.c.) , where 
A v is the superconducting order parameter (An = 0). The tunnelling Hamiltonians are 
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Htunn,r) = J2ka (^fer^cr + H.c.). Here, V v are the spin- and wavevector-independent 
tunnel matrix elements, and c^^cL.^) and d a (d\) represent the annihilation (creation) 
operators for the leads and dot, respectively. The tunnel-coupling strengths are 
characterized by T v = 2n\V ri \ 2 J2k ~ e fc)- 



3. Current formula 



We start with deriving a general formula for the charge current in lead 77 by using 
the approach of Ref. [29] generalized to superconducting leads. Similar formulae 
that relate the charge current to the Green's function of the dot in the presence of 
superconducting leads have been derived in previous works, in particular for equilibrium 
situations, see e.g. Refs. [TO], [22j. The formula derived below is quite general, as 
it allows for arbitrary bias and gate voltages, temperatures, and superconducting 
order parameters for a quantum dot coupled to an arbitrary number of normal and 
superconducting leads. For this, it is convenient to use the operators ip v k = (c v hi, ct fc|) T 
and 4> = (d-f,c?t) T in Nambu formalism. The current from lead n is expressed as 
J v = e(dNjdt) = i(e/h)([H,N v }) = i(e/h)([H tmn>v , N v }) Q, with N v = Efc^Jfc^fc, 
where T\ , Ti , r 3 indicate the Pauli matrices in Nambu space and e > the electron 
charge. Evaluating the commutator leads to 

J n = -f E/ ^Re{Tr [r 3 V r? G< D ,, fc M]}, 

k ^ 

r rj = Diag(K ? , — V*) and the lead-dot lesser Green's functions (G < r>, v k(^)) J 



(1) 



with = Diag(l^, —V*) and the lead-dot lesser Green's functions (G < i 
that are the Fourier transforms of «(^fc n (O)0 m (t)). In the following, we assume the 
tunnelling matrix elements V v to be real (any phase of V v can be gauged away by 
substituting A,, — > A,, exp(— 2i arg V v )). The Green's function G < D, J7 fc is related to 
the full dot Green's functions and the lead Green's functions by a Dyson equation in 



Keldysh formalism: G< D , vk {uj) = G R (cu)Vtg< (u) + G< 
the retarded (lesser) dot Green's function, and and g 



A (<), 

r)k 1 



MVtg^(u; 



where G 11 ^ (to) is 
the lead advanced (lesser) 



Green's function. Using this relation and assuming energy-independent tunnel rates T^, 



we obtain for the current J„ 



JlT) + ^2r? With 



J 



1'/ 



du 

—T n D v (uj)lm 
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[i- A ' I 









(2G»/» + g<m)]} 
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2,, 



|/^r^(o;)Re 



Tr 



A, 



G< 



10J 



(2) 
(3) 



where A r 



and fr,{y>) — [1 + exp(c<j — ji rj ) j '(k^T)] 1 is the Fermi 



' 
I a; o . 

function, with T being the temperature and fee the Boltzmann constant. The dot 
Green's functions (G < i)(uj)) m _ and ( G R d(i^) ) are defined as the Fourier transforms 

''m,n V / m,n 



% Note that [H V> N V ] + but ([H v> N v ]) = for 77 = S l ,r. 
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of i{(fil l (0)(j) m (t)) and — i9(t)({(f) m (t), 0^(0)}), respectively. The two weighting functions 
D v {uj) and S v (uj) are given by 

D v (u) = 0(\u - N \ - |A„|) 



u - - |A 



If \ 



S v (uj) = , ^ Q(\W - \oj- Hr,\)- 



\A V \ 2 -(u- fin) 2 

The terms Ji v and involve excitation energies ui above and below the superconducting 
gap, respectively. For 77 = N, only the part of J\ v that involves normal (diagonal) 
components of the Green's functions contributes, and the current reduces to the result 
presented in [29] . For superconducting leads, this part describes quasiparticle transport 
that is independent of the superconducting phase difference. The other part of J\ v 
involves anomalous (off-diagonal) components of the Green's functions and is, in general, 
phase dependent. The contribution to the Josephson current stemming from this term is 
the dominant one in the regime considered in [T3J Q3J [15] . The excitation energies above 
the gap are only accessible either for transport voltages exceeding the gap or by including 
higher-order tunnelling, involving virtual states with quasiparticles in the leads, and, 
therefore, J\ n vanishes for large \A V \. In this case J2 V , that involves only anomalous 
Green's functions with excitation energies below the gap, dominates transport. It is, in 
general, phase dependent, and describes both Josephson as well as Andreev tunnelling. 
In the following we consider the limit \A V \ — > 00, where the current simplifies to 

Jr, = ^r, J Kd ; d T )|sin(*-^), (4) 

with being the phase of A,, and = |(d|df)| exp(i\l/) the pair amplitude of the 

dot that has to be determined in the presence of Coulomb interaction, coupling to all 
(normal and superconducting) leads and in non-equilibrium due to finite bias voltage. 

We now consider a symmetric three-terminal setup with Ts L = Ts R = Tg, 
A Sl = |A|exp(z$/2) and A Sr = |A| exp(— z$/2), and = /is R = 0. The quantities 
of interest are the the current that flows between the two superconductors (Josephson 
current) Jj OS = (Js L ~ ^s R )/2 and the current in the normal lead (Andreev current) 

Jand — Jn = ~(Js L + ^S R )- 

Furthermore, we focus on the limit of weak tunnel coupling, Ts < k^T. In this 
regime, an Josephson current through the dot in equilibrium would be suppressed even 
in the absence of Coulomb interaction, U = 0, since the influence of the superconductors 
on the quantum-dot spectrum could not be resolved for the resonance condition e f» 0. 
This can, e.g., be seen in the exactly-solvable limit of U = together with T N = 0, 
where the Josephson current is Jj OS = (e/2ft)rf sin($) [/(-e A ($)) - f(ex(^))] Aa($) 



with e A ($) = ye 2 + T| cos 2 ($/2). This provides an additional motivation to look for a 
non- equilibrium mechanism to proximize the quantum dot. 



N on- Equilibrium Josephson and Andreev Current through Interacting Quantum Dots 6 



4. Kinetic equations for quantum-dot degrees of freedom 

The Hilbert space of the dot is four dimensional: the dot can be empty, singly 
occupied with spin up or down, or doubly occupied, denoted by |%) G {|0), | t), 
| i),\D) = d\d\\0)}, with energies Eq, E^ = E\, Er>. For convenience we define the 
detuning as 5 = Ed — E = 2e + U. The dot dynamics is fully described by its reduced 
density matrix p D , with matrix elements P* 1 = (pd)x2xi- The dot P a ^ r amplitude (e^df) 
is given by the off-diagonal matrix element Pp. The time evolution of the reduced 
density matrix is described by the kinetic equations 

jPtm + {{E X1 - E X2 )P£(t) = £ t dt'W^(t,t')P^(t'). (5) 

xi,x' 2 

We define the generalized transition rates by W*™£ = J-oo ^'^x^xf wn i° n are tne 
only quantities to be evaluated in the stationary limit. Together with the normalization 
condition S X P X — 1, (151) determines the matrix elements of pd- Furthermore, in ([5]) 
we retain only linear terms in the tunnel strengths T v and the detuning S. Hence, we 
calculate the rates W*^ 1 to the lowest (first) order in T v for 5 = 0. This is justified in 
the transport regime T s , T N , 5 < k B T. 

The rates are evaluated by means of a real-time diagrammatic technique |30j . 
that we generalize to include superconducting leads. This technique provides a 
convenient tool to perform a systematic perturbation expansion of the transport 
properties in powers of the tunnel-coupling strength. In the following, we concentrate 
on transport processes to first order in tunnelling (a generalization to higher orders is 
straightforward). This includes the transfer of charges through the tunnelling barriers 
as well as energy-renormalization terms that give rise to nontrivial dynamics of the 
quantum-dot degrees of freedom. 

We find for the (first-order) diagonal rates W XlX2 = the expressions W a Q = 

r N / N (-C//2); W a = r N [l - f N (-U/2)];W Da = T N f N (U/2);W aD = T N [l - fa(U/2)]. 
The N lead also contributes to the rates W° D = {W° D ° D )* = -T N 1 + f N (-U/2) - 
h{U/2) + zB] where B = ±Re [ij> (§ + z^gf ) - ip (f + , with // N being 

the chemical potential of the normal lead and if)(z) the Digamma function. Notice that 
B vanishes when //^ = or U = 0. The superconducting leads do not enter here 
due to the gap in the quasi-particle density of states. These leads, though, contribute 
to the off-diagonal rates W °g = Wg? = (W^)* = (w%)* = -W^° = -W£ D = 

- (w™)* = - (wS D Y = -*r s cos($/2). 

For an intuitive representation of the system dynamics we define, in analogy to |31j . 
a dot isospin by 



_ P D + P° D . _, P D -P^ . . _ Pp - Po 



(6) 
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From (151), we find that in the stationary limit the isospin dynamics can be separated 
into three parts, = dl/dt = (dl/dt) &cc + (dl/dt) icl + (dl/dt) rot , with 

f ) = - ^[1 - M-U/2) - h(U/2)]e z (7) 
= - r N [l + M-U/2) - / N (C//2)]I (8) 
= I x B eff (9) 

rot 

where e z is the z-direction and B e fj = {2Fs cos($/2), 0, — T^B — 2e — U} is an effective 
magnetic field in the isospin space. The accumulation term (jTJ) builds up a finite isospin, 
while the relaxation term (jSJ) decreases it. Finally, ([9]) describes a rotation of the isospin 
direction. 




5. Non-equilibrium Josephson current 

In the isospin language the current in the superconducting leads is 

Js L , R = ^r s [I y cos($/2) ± I x sin($/2)] , (10) 

where the upper (lower) sign refers to the left (right) superconducting lead. The I y 
component contributes to the Andreev current, while I x is responsible for the Josephson 
current. To obtain subgap transport, we first need to build up a finite isospin component 
along the z-direction, i.e. we need a population imbalance between the empty and doubly 
occupied dot [(this is generated by the accumulation term in (J7J)]; second, we need a 
finite B e fj which rotates the isospin so that it acquires an inplane component. In order 
to have a finite Josephson current (I x ^ 0), we need the z-component, — T^B — 2e — U, 
of the effective magnetic field producing the rotation to be non zero. 
The Josephson current and the Andreev current read 
er s [2e + f/ + r N 5]r s sin(<l>) 



j, 



and 



jos n |B crT | 2 + r2 J [i + / N (-f//2)-/ N (f//2)] 2 

l-f N (-U/2)-f N (U/2) 
1 + M-U/2) -h(U/2) 
e£s 2r N r s [l + cos(<I>)] 

n |B eff |2 + r 2 <[i + / N (-[//2)-/ N ([//2)p 

x [1- M-U/2) -MU/2)). (12) 

These results take into account only first-order tunnel processes, i.e. the rates 
are computed to first order in T^. The factor [1 — f N (—U/2) — f N (U/2)} ensures that 
no finite dot-pair amplitude can be established if the chemical potential of the normal 
lead, /xn, is inside the interval [— U/2, U/2] by at least k^T. In this situation both 
the Josephson and the Andreev currents vanish. On the other hand, this factor takes 
the value — 1 if /jn > U/2 and the value +1 if /jn < —U/2. Hence, the sign of the 
Josephson current can be reversed by the applied voltage (voltage driven 7r-transition). 
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The considerations above establish the importance of the non-equilibrium voltage to 
induce and control proximity effect in the interacting quantum dot. In figure [2] we show 
in a density plot (a) Jj OS and (b) J an( j for $ = n/2 as a function of the voltage /in and 
the level position e. Both the control of proximity effect by the chemical potential /iN 
and the voltage driven 7r-transition are clearly visible. If the detuning is too large, 
1 5 + T N B\ > + 4r§ cos 2 ($/2), it becomes difficult to build a superposition of 

the states |0) and \D), which is necessary to establish proximity. As a consequence, 
the Josephson and the Andreev current are algebraically suppressed by <5 _1 and 5~ 2 , 
respectively. Figure [3] shows the Josephson current as a function of 5 = 2e + U. The 
fact that the Josephson current is non zero for 5 = is due to the term T-^B, i.e. of the 
interaction induced contribution to the ^-component of the effective field B e g- acting on 
the isospin. The term \B\ has a maximum at /xn = U/2, which causes this effect to be 
more pronounced at the onset of transport. The fact that the value of the Josephson 
current varies on a scale smaller than temperature indicates its nonequilibrium nature. 

A 7r-transition of the Josephson current can also be achieved by changing the sign 
of 5 + r N I?, as shown in figure H] where Jj OS is plotted as a function of the phase difference 
$ for different values of the level position. Notice that the current for 5 = (e = —U/2) 
is different from zero only due to the presence of the term T^B acting on the isospin. 

6. Conclusion 

In conclusion, we have studied non-equilibrium proximity effect in an interacting 
single-level quantum dot weakly coupled to two superconducting and one normal lead. 
We propose a new mechanism for a Josephson coupling between the leads that is 
qualitatively different from earlier proposals based on higher-order tunnelling processes 
via virtual states. Our proposal relies on generating a finite non- equilibrium pair 
amplitude on the dot by applying a bias voltage between normal and superconducting 
leads. The charging energy of the quantum dot defines a threshold bias voltage above 
which the non-equilibrium proximity effect allows for a Josephson current carried 
by first-order tunnelling processes, that is not suppressed in the limit of a large 
superconducting gap. Both the magnitude and the sign of the Josephson current 
are sensitive to the energy difference between empty and doubly-occupied dot. A tt- 
transition can be driven by either bias or gate voltage. In addition to defining a threshold 
bias voltage, the charging energy induces many-body correlations that affect the dot's 
pair amplitude, visible in a bias-voltage- dependent shift of the 7r-transition as a function 
of the gate voltage. 
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Figure 2. Density plot of the a) Josephson and b) Andreev current, for fixed 
superconducting-phase difference $ = tt/2, as a function of the dot-level position e 
and of the chemical potential of the normal lead /iN- The symbols ± refer to the sign 
of the current. The other parameters are T s = T N = 0.01C/, and k B T — 0.05C7. 
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Figure 3. Josephson current, for fixed superconducting-phase difference <£> = ir/2, as 
a function of the detuning 5 — Ed — Eq = 2e + U for different values of the chemical 
potential. The other parameters are = Tn — 0.01J7 and k^T — 0.05J7. 




Figure 4. Josephson current as a function of the superconducting-phase difference 
$ for different values of the detuning. The other parameters are r$ = Tn = 0.01C/, 
/i N = U, and k B T = 0.05U. 
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